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We investigate a potentially large CP violating asymmetry in the neutral Higgs 
boson decay into a heavy quark pair or a W + -W~ pair. The source of the CP 
nonconservation is in the Yukawa couplings of the Higgs boson which can contain both 
scalar and pseudoscalar pieces. One of the interesting consequence is the different 
rates of the Higgs boson decays into CP conjugate polarized states. The required 
final state interactions can be either the strong or the electroweak interactions. The 
CP violating asymmetry can manifests itself in the asymmetry in the energies of the 
secondary decay products of these heavy quarks and W bosons. Such asymmetry 
can be measureable in the future colliders such as SSC or LHC. 



PACS numbers: 11.30.Er, 12.15.Cc, 14.80. Gt 



Typeset Using REVTEX 



The physics related to the Higgs boson is the most mysterious part of the Standard 
Model. It is widely expected that whatever information we can obtain about the Higgs 
sector from the next generation colliders will give us hints about the potential new physics 
beyond the Standard Model. Most of analyses about the search for Higgs bosons, in the 
Standard Model or beyond, often ignore potential CP violation. However, in these models, 
the lightest neutral Higgs boson can have interesting CP violating phenomena. In some 
sense, as we shall elaborate later, the CP violating aspect may be the most interesting part 
of the Higgs boson physics beyond the Standard Model once a neutral Higgs is identified. 

In this letter, we investigate an interesting signature of CP violation in the neutral Higgs 
decay into a heavy quark pair QQ or a W + W~ pair. For the heavy quark QQ mode, the 
polarizations of the quarks are either QlQl or QrQr. (Note that we have adopted the 
notation that Ql is the antiparticle of Q R and should be left handed.) These two modes 
happen to be CP conjugate of each other. Therefore one can immediately consider the CP 
asymmetry in the event rate difference, N(Q L Q L ) — N{Q R Q R ). For the W pair mode, the 
two polarization modes W^W^ and H^W^r are also CP conjugate of each other. (Here 
we denote the states of helicities —1, 1 and by L, R and || respectively.) Therefore, their 
difference in the decay rates, N(W^W£) — N(W R W R ), signals CP violation. For the Z 
pair mode, one can look at the similar asymmetry N(ZlZl) — N(Z R Z R ) also. 

A method in detecting the asymmetry N(Q L Q L ) — N{Q R Q R ) was proposed recently in 
Ref . . One assumes that the quark Q decays semileptonically through the usual V — A 
weak interaction. For a lighter quark (such as the b quark), the information about polar- 
ization is probably washed away in the soft process of hadronization before its subsequent 
decay. However, for a heavy quark such as the top quark, since the hadronization time is 
much longer than the decay time ||, one can analyze polarization dependence of its de- 
cay at the quark level. The top quark first decays into a b quark and a W + boson, which 
subsequently becomes Iv. For heavy top quark, the W + boson produced in top decay is 
predominantly longitudinal. Due to the V — A interaction, the b quark is preferentially pro- 
duced in the left-handed helicity. So the longitudinal W + boson is preferentially produced 

2 



along the direction of the top quark polarization. Therefore the anti-lepton I produced in 
the W + decay is also preferentially in that direction. At the rest frame of t, the angular 
distribution |3J of the produced I has the form 1 + cosip, with ip as the angle between I and 
the helicity axis of t. When the Higgs boson decays, the top quark is produced usually with 
nonzero momentum. As a result of the Lorentz boost, the anti-lepton I produced in the 
decay of the right handed top quark t R has a higher energy than that produced in the decay 
of the left handed top quark tt- Similarly, the I lepton produced in the decay of it has a 
higher energy than that produced in the decay of Ir. Consequently, in the decay of the pair 
tjJt the lepton from it has a higher energy than the anti-lepton from tt] while in the decay 
of t^tn the anti-lepton has a higher energy. Therefore one can observe N(ttit) — N(t R in) 
by measuring the energy asymmetry in the resulting leptons ||. 

For asymmetry in the W + W~ mode, one can look at the leptonic decays of the W gauge 
bosons. We are interested in the transverse W bosons in this case. In the rest frame of W + , 
which decays into Iv, the angular distribution of I has the form (1 + costp) 2 , with ip as the 
angle between / and the helicity axis along which the spin projection of W is one. Similar 
to previous analysis, the anti-lepton in the decay of W R has a higher energy than the lepton 
from the Wr, while the lepton in the decay of W£ has a higher energy than the anti-lepton 
from the W£. Therefore, the relative energy between the lepton and the antilepton in the 
decay of a W + W~ pair can be translated into the information about the polarization of W 
bosons. The asymmetry is much harder to detect in the leptonic final states for the case of 
Z pairs. Because of the dominance of the axial-vector coupling to leptons, there is only a 
weak angular dependence for the outgoing leptons. However it is present in principle. 

In order to generate the asymmetries N(t L i L ) - N(t R i R ) or N(W^W£) — N (W R W R ) , it 
is necessary to include effects of the final state interactions in order to escape from the her- 
miticity constraint at the tree level due to the CPT theorem. These final state interactions, 
as shown in Fig. 1, can come from the strong, the electroweak, or the Higgs corrections to 
the HQQ or HW + W~ vertices. 

In the H — > ti channel, CP non-conservation occurs in the complex Yukawa couplings, 
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C-mt = -(m t /v)t{AP L + A* P R )tH , v = {V2G F )~* ~ 246 GeV . (1) 

The complex coefficient A is a combination of model-dependent mixing angles. Simultaneous 
presence of both the real part A R = KeA and the imaginary part Aj = 1mA guarantees CP 
asymmetry. For example, at the low energy regime, it can give rise to the electric dipole 
moment of elementary particles |||7| . Here we will show that CP nonconservation manifests 
itself in the event rate difference in collider experiments. We denote 5^ CD , 5 1 , 5 Z , 8 , 
8 WW ; 8 ZZ , as contributions in CP asymmetry due to the exchange of the gluon, the photon, 
the Z boson, the Higgs bosons, or those with intermediate states of a W pair or a Z pair 
respectively, 

N(t L t L ) - N(t R i R ) = £Qcd + 6 y + s z + 6 H + 5 WW + 5 ZZ 
N(al\ tt from H) f3?A 2 R + Aj ' U 

with 01 = 1 — Amf/Mfj. The one-gluon exchange gives large CP asymmetry, 

8^ CD = (Ca s /2)Im(A 2 )(l-i3^, (3) 

with the color factor C — 4/3. It is interesting to note that the imaginary part of the 
one-loop graph contributes a factor of A R while the tree graph contributes Aj (that is, the 
pseudoscalar coupling). When Aj ~ A R , the asymmetry is of order of the strong coupling 
as, about 10 _1 . Fig. 2 shows how such asymmetry depends on Mh- Note that there is no 
strong constraints ||,[7| on lm(A 2 ), which can easily be of order one. 

The electromagnetic correction <5 7 is obtained by replacing as by the QED coupling a, 
and C by 4/9, the charge squared of the t-quark. The result is negligible. The contribution 
by the Z-exchange is 



8 



z 



^Mfl{P(f) - W + g )H + (1 - W3«} . (4) 



Here the function F(x) = 1 — x _1 log(l + x), and x\y = sin 2 9\y ~ 0.23. The term with 
the factor (1 — 8xw/3) 2 is due to the part of the one-loop graph in which both vertices 
of Z boson are vectorial while the rest are due to the part in which both vertices of Z 
are axial-vectorial. The parts with one axial vector vertex and one vector vertex do not 



contribute. For the part involves only the axial vector vertices of Z boson, the tree graph 
can contributes either scalar or pseudoscalar coupling. Contribution from the Higgs bosons 
is more complicated. In general, there are more than one neutral Higgs bosons ifj, with 
Yukawa couplings to the t-quark of the form, C Hi tt — —(jfk/v) Y,ii(a-iPL + a*Pn)tHi. The 
coefficient A is only one of these a,. The overall effect from all Higgs bosons is 

S H = E Fi^ 1 ) [$MA 2 *?) + (1 - A 2 )Im(A)Re(AH 2 + Aaf)] . (5) 

i Hi 

This general formula is valid for the decay of each Higgs boson in the system. For the 
purpose of illustration, we work in a simplified case of the decay of the lightest Higgs boson 
H. Then, if the effects of heavier Higgs bosons can be neglected, the formula is much reduced 
as an = A and the first term in the bracket of the above equation vanishes. Numerical study 
shows that 5 H is very small in this scenario. 

There are important contributions from processes involving W + W~ or ZZ intermediate 
states as shown in Fig. 1. The amplitudes depend on the the vertices, 

Bg 2 H[M w W +v W~ + \{M z /cos9 w )Z v Z v \ . (6) 

In any model, the scalar boson that couples to the W pair or Z pair is the scalar partner of the 
unphysical Higgs boson. This scalar boson is in general not a mass eigenstate. Therefore we 
parametrize its coupling by a phenomenolgical coefficient B in addition to the usual SUl{2) 
gauge coupling g 2 . The coefficient B represents the mixing factor needed to reach the mass 
eigenstate. It reduces to unity in the Standard Model with one single Higgs doublet. For 
the diagrams with W + W~ intermediate states, we obtain the CP asymmetry, 



. ww _ ^GpM^BAjPwPt r r 3& 2 -20*y-l , 1 - ft 
An ' 

-2L(^,A)+2 + i±|f}, (7) 
1 ~ Pw 



f r 3ff-2/fr-l 1-ft n 



where (3^ = 1 — AM^/Mjj, and the function 



x 2 — y 2 



x-y 
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Also, for the diagrams with ZZ intermediate states, we obtain the CP asymmetry, 
r „ _ V2G F MIBA^ Z r _ 8 + 1±^ + M J0(1 _ ^ 



16tt# L v 3 v 4^ 



i-A 2 (i + /?l) 2 , i + ^I 



with 



K = lot 



_i ' t V ' i /j / jy- _ i _ q ' 

1 2/9zA 1-/3! 



1 + /?1 + 2A/3 Z 



(10) 



'z 

The first term in the square bracket is the contribution where both of the Z couplings with 
the top quark are vectorial. The remaining two terms correspond to the contribution where 
both couplings are axial- vectorial. Also note that the one-loop induced vertex is always 
scalar and the Aj factor arises from the tree level amplitude. Numerical study shows that 
fiWW an( j fizz are ver y important for the natural scenario Ar ~ B. Their contributions 
can dominate when Mh is large as illustrated in Fig. 2. In this letter, we shall ignore the 
contributions due to tri-Higgs couplings because they are more model dependent in general. 
Also, this type of contributions disappears, due to the vanishing imaginary part, when the 
decaying Higgs boson is the lightest Higgs boson. In any case, this type of contribution will 
not affect our results very much as long as no accidental cancellation occurs. 

Since CP violation in Eq.(l) would disappear if the top quark were massless, one may 
wonder why the above formulas for 5 ZZ and S ww are not proportional to m t . The answer 
is that, in the denomenator, the leading (tree) level contribution also has a factor of m t is 
our parametrization in Eq.(l). The m t required physically in the numerator is artificially 
cancelled by this factor in the denomenator. 

The polarization asymmetry is Eq. (0) can be translated in to the lepton energy asymme- 
try |l];|2|f|. The energy Eq(1 + ) distribution of a static t quark decay t — > l + ub is very simple 
Q in the narrow width T\v approximation when is negligible. 

x (l-x )/D if m 2 w /m 2 t < x < 1, 
f(xo) = { (11) 
otherwise. 
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Here we denote the scaling variable x$ = 2Eo(l + )/m t and the normalization factor D = 
| — \ {mw /m t ) 4 + |(miy/m 4 ) 6 . When the t quark is not static, but moves at a speed (3 with 
helicity L or R, the distribution expression becomes a convolution, 

fRM x iP)= f( x o) — Trrm — dx o ■ ( 12 ) 

Here x = 2E(l + )/E t . The kernel above is related to the (1 ± cosip) distribution mentioned 
in the introduction. Similar distributions for the t decay is related by CP conjugation at the 
tree-level. Using the polarization asymmetry formula in Eq.(||), we can derive expressions 
for the energy distributions of l~ and l + : 

N-'dN/dx^) = ~(1 ± A)f L (x,(3 t ) + |(1 T A)/ fl (z, A) • (13) 

Here distributions are compared at the same energy for the lepton and the anti-lepton at 
the H rest frame, x(l~) = x(l + ) = x = 4£'(/ ± )/M^. To prepare a large sample for analysis, 
we only require that each event has at least one prompt anti-lepton Z + from the t decay or 
one prompt lepton l~ from the t decay. Fig. 3 compares the lepton and anti-lepton energy 
distribtuions for a typical asymmetry A = 0.1. One can also sum over channels of different 
lepton flavor to increase the event rate. 

For the W pair production, the tree level amplitude, parametrized by B, interferes with 
the fermion-loop amplitude (Fig. 4) of the b quark exchange to produce the CP violating 
asymmetry. The resulting asymmetry is 

_ N{WjWE) - N(W£Wn) _ 3y/2GW ft , 2(1 - (3 2 w ) 2 \ 

Aw = iV(all W+W-, from H) " AnB Pw ^3 - 2^ w + 3(3^ ' (14) 

Note that, for a heavy Higgs boson, its leading mode is neither W^W£ nor Wr W^, but the 
longitudinally polarized state of . However, the W^Wm - state is CP self conjugate and 

cannot provide the CP asymmetry we are looking for. Consequently, there is a suppression 
factor in the last bracket of Eq. (|I4|) . The factor is the ratio between the transverse W~W + 
rate and the overall W~W + rate from the Higgs boson decay. Fig. 5 shows such asymmetry, 
which is sizable near the threshold of ti. As expected, when the ti threshold gets higher the 



asymmetry gets smaller due to the suppression of the transverse modes mentioned above. 
This asymmetry can be detected by measuring the energy asymmetry of the oppositely 
charged leptons in the decays of the bosons. We define the dimensionless variable 
x(l) = 4E(1)/Mh in the H rest frame as before. The x distributions due to the process 
H — > W + W~ are given by 

1 dN (1 + P W ) 2 n3[^-(1-i) 2 ] 



Ndx{l±) V 3 _ 2/3^ + 3{3tv J 4/3, 3 



w 

2 \2 



v / (1-(3 W ) 2 sA w ^ 3(x-l±s(3 w ) 2 
s= k+i [ 3 - + Wtv 2 ) 8^ 1 " J 

The first term above comes from the longitudinal W, and the rest are due to the two trans- 
versely polarized W's (helicity s = —1,-1-1). Each bracket corresponds to the production 
weight of the corresponding polarized state. Fig. 6 compares the lepton and anti-lepton 
energy distribtuions for a typical asymmetry Ay/ = 0.01. As before, we only require that 
each event has at least one prompt anti-lepton l + from the W + decay or one prompt lepton 
l~ from the W~ decay. It is worthwhile to point out that the smallness of Aw is mostly due 
to the large contribution of the longitudinal W mode. For example, this suppression factor, 
in the last bracket of Eq.(0), is 0.015 when M H = 400 GeV. In Fig. 6, We expect that 
events with the lepton energy x ~ (l±/3w) come from the transverse modes, and otherwise, 
the central bump is dominated by the longitudinal W mode. Therefore, in principle, one can 
enhance effects from the transverse modes by imposing energy cuts that eliminate events in 
the central region from the longitudinal mode. The efficiency of such cuts deserves detailed 
simulations in the future p. 

Similar CP asymmetry can be formed between the ZlZl and the ZrZr pairs, 



N(Z L Z L ) - N(Z R Z R ) = 3V2G F m 2 t ft r K 8 2 Kf3 z 

N(Z L Z L ) + N(Z R Z R ) AnB W Wz 1 3 XW) *Pt 



However it is not clear how feasible one can decode such asymmetry from the final decay 
products of the Z's. 

In conclusion, we have shown a potentially large CP violating asymmetry in the decay of 
the neutral Higgs boson. For the channels H — > ti, W + W~, ZZ, to be available, the Higgs 
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mass has to be about 180 GeV or higher. That rules out LEP II and possibly Tevatron as the 
production machine. The heavy Higgs boson can be produced copiously in the future hadron 
collider via the gluon-gluon or W-W fusions, or in the next high energy e + e~ machine via 
Z-bremsstrahlung or W-W fusion. At SSC, if 200GeV < M# < ITeV, there will be between 
10 6 and 10 5 Higgs boson produced |TD[ for an annual integrated luminsity of 10 4 pb _1 . The 
"gold-plated" events, H — > ZZ — > (/ + Z~)(/ + Z _ ), are supposed to be the discovery mode 
for the Higgs mass up to about 600 GeV. The W + W~ mode is about twice as prolific as 
the ZZ mode and has even larger leptonic branching fraction. However, their signature of a 
charged lepton plus two jets from the W has large non-resonant background in the Standard 
Model. More works || are needed in decoding the signature. The branching fraction for 
the Higgs boson decay into tt can be between 5% and 20% for 350GeV < M H < ITeV and 
m t = 150GeV for example. Their signature is just as hard to decode as that of the W + W~ 
mode. However, once the Higgs mass is known from the ZZ mode, the W + W~ and tt modes 
should be much easier to identify. 
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FIGURES 

FIG. 1. Feynman diagrams for the process H — > ti. The tree amplitude (a) interferes with 
those one-loop amplitudes with the final state interactions coming from (b) the exhange of a gluon 
g, a photon 7, a Z gauge boson or a Higgs boson Hi; (c) the intermediate W + W~ boson pair; (d) 
the intermediate ZZ boson pair. 

FIG. 2. [N(t L i L )-N(t R i R )]/N(al\ ti, from H) versus M H for m t = 150 GeV. The illustration 
is for the case Aj = A R = B = 0.5. 

FIG. 3. The energy distribution of the anti-lepton (lepton) in the H rest frame is plotted in 
the solid (dashed) curve, for the case that H -> it, M H = 300 GeV, m t = 120 GeV, and A = 0.1. 
The asymmetry is shown in the bottom. 

FIG. 4. (a) The one-loop diagram for H — > W + W~ produces the CP violating asymmetry 
when it interferes with the tree amplitude, (b) The same effect occurs for H — > ZZ. 

FIG. 5. [N(W£W£) - N(W£\V R )]/N(W + W- , from H) versus M H for different values of 
m t =100 GeV (solid), 150 GeV (dashed), 200 GeV (dotted), and 250 GeV (dash-dotted). 

FIG. 6. Energy distributions of the lepton and the anti-lepton in the H rest frame are 
compared, for the case that H — > W + W~, Mh = 300 GeV, and A\y = 0.01. The asymmetry is 
shown in the bottom. 
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